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We investigate cross-relaxation interactions between Tm and A1 in Tm 3+ :YAG using two optical 
methods: spectral holeburning and stimulated echoes. These interactions lead to a substantial 
reduction in the hyperhne lifetime at magnetic fields that bring the Tm hyperfine transition into 
resonance with an A1 transition. We develop models for the measured echo decay curves and 
holeburning spectra near a resonance, which are used to show that the Tm-Al interaction has a 
resonance width of 12 ± 3 kHz and reduces the hyperfine lifetime at resonance to 0.5 ± 0.3 ms. 
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I. INTRODUCTION 

Rare earth ions in crystals have the unique property 
of maintaining long coherence times at high spatial den¬ 
sity. This property lies at the root of the sustained in¬ 
terest in these materials for information processing ap¬ 
plications in both the classical and quantum domains. 
When combined with large inhomogeneous linewidths, 
long coherence times allow classical processing with very 
large time-to-bandwidth products, and applications of 
rare earth materials to RADAR, electronic warfare, and 
astrophysics have been pursued actively for the last two 
decades (Ref Q] and references therein). These applica¬ 
tions favor materials such as Tm 3+ :YAG (Y3AI5O12) and 
Er 3 +:Y 2 Si0 5 as they have a one-to-one correspondence 
between the signal frequency and the storage ion transi¬ 
tion frequency. 

For quantum applications, the long coherence times 
of the easily accessed rare earth optical transitions al¬ 
ready make them good candidates, but the possibility of 
transferring optical coherence to much longer lived spin 
coherence makes these materials particularly attractive. 
The non-Kramers ions Pr, Eu, and Tm in low symmetry 
sites, such as in YAG and Y 2 SiOs, are particularly useful, 
as they have singlet electronic levels that are highly in¬ 
sensitive to the crystalline environment, a property that 
favors the existence of very long coherence times. The 
quantum applications of rare earth ions are reviewed in 
Refs El and G2 

Optimising performance for these applications requires 
a detailed understanding of the way the rare earth ion 
interacts with its crystalline environment. Ion-ion in¬ 
teractions affect properties such as the optical density, 
inhomogeneous linewidth, and level lifetimes, and are of¬ 
ten the dominant source of decoherence. For example, in 
Er 3 +:Y 2 Si05, Bottger et al. showed that interactions be¬ 
tween Er ions dominate optical dephasing and that this 
effect can be avoided by using low concentration crys¬ 
tals in high magnetic fields, resulting in a coherence time 
of 2 ms, the longest optical coherence time observed in 
a solid 4 . For hyperhne transitions of non-Kramers ions, 


meanwhile, dephasing usually comes from the magnetic 
fluctuations due to hipping lattice spins, and suppressing 
this effect contributed to the longest hyperhne coherence 
times measured, of 6 hourtP. 

Despite those unquestionable successes, the dynamics 
of the host crystal, a dominant source of decoherence, 
have still to be better understood. In the present work we 
intend to gain information on the dynamics of lattice ions 
through their resonant interaction with optically active, 
low-concentration, impurities. In this paper, we use the 
term “resonant interaction” to refer specifically to spin 
polarization transfer between the interacting ions. 

Heteronuclear resonant interactions, arising from the 
magnetic dipole-dipole coupling of unlike nuclei, are well 
known in NMR spectroscopy. They are used in a vari¬ 
ety of common techniques, including nuclear Overhauser 
spectroscope®, cross polarisation spectroscopy®, and adi¬ 
abatic demagnetisatioiP, for detecting low concentration 
or low magnetic moment species by transferring spin 
polarization from a high concentration, high moment 
species. The situation we consider is just the opposite. 
The high sensitivity of optical detection enables us to 
directly monitor scarce impurity ions. The lattice dy¬ 
namics are little affected by those ions, but can be seen 
in their time evolution. We take advantage of resonance 
to precisely select the lattice ions that are probed by the 
impurities. 

Techniques based on resonant interactions between un¬ 
like nuclei have very seldom been used in optical spec¬ 
troscopy, particularly in the high resolution spectroscopy 
used to study rare earth materials. We are only aware 
previous work on one material, P^LaF^EM- 1 ®. In this ma¬ 
terial, strong resonant interactions between Pr and F oc¬ 
cur. These were used to prove the existence of a frozen 
core of fluorine spins®. 

In this paper we investigate resonant interactions be¬ 
tween Tm and its Al neighbors in Tm 3+ :YAG, a material 
considered for both cl assi cal and quantum information 
processing applications 1 These interactions cause 

sharp decreases in the hyperfine state lifetime when a 
magnetic field brings the two ions into resonance. This 
occurs for fields < 10 mT, a region commonly used for 
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information processing applications. High temperature 
NMR investigation of the resonant exchange of excitation 
between Tm and A1 in the isomorphic crystal TmAIG 
has been undertaken in the pasiP^. It was found that the 
resonant interaction with the abundant Tm ions can fa¬ 
cilitate spin-lattice relaxation of 2 ' Al. In this paper, we 
work with a low concentration of Tm, precisely so that 
the Al dynamics are not altered. 

Our investigation is based on two common optical spec¬ 
troscopy methods, spectral holeburning and stimulated 
echoes. We start by giving examples of how ion-ion reso¬ 
nant interactions affect the holeburning spectrum in Sec¬ 
tion HU In Section |III[ we characterize the relaxation be¬ 
havior of the system away from resonance using stimu¬ 
lated echoes. We then use this information to develop a 
simple model for the echo decay at resonance. In Sec¬ 
tion [TVAj we extend this model to the more complicated 
case of the holeburning structure. In Section |IVB[ we 
include the Al Zeeman interaction into this model, and 
in Section [IV C[ compare the model to the experimental 
holeburning data. Finally, in Section[V]we briefly discuss 
the resonant interaction mechanism, and the possibility 
of directly modelling the dynamic transfer of the Tm ex¬ 
citation to the spin bath. 


II. HOLE BURNING EXTINCTION DUE TO 
RESONANT INTERACTIONS 

Tm has a single stable isotope, 169 Tm, with a nuclear 
spin of I = In Tm 3+ :YAG, the Tm ion occupies the 
single Y site of D 2 symmetry. It has an optical transition 
from the 3 H 6 (0) ground state to the 3 H4(0) excited state 
at 793 nm. In both the optical ground and excited states, 
the mixing of the electronic wavefunctions by the crystal 
field leads to a large and extremely anisotropic enhanced 
nuclear Zeeman tensoiG^E^of the order of 100 MHz/T. 

Resonant ion-ion interactions have a particularly 
striking effect on the holeburning spectra of the 
3 Hg(0) —)- 3 H4(0) transition of Tm 3+ :YAG. In zero mag¬ 
netic field, Tm 3+ :YAG has a holeburning lifetime of 
~ 10 ms due to storage in the 3 F 4 metastable state. 
The lifetime increases to several minutesP-^l when a mag¬ 
netic field is applied as this splits the doubly degenerate 
ground state, allowing storage in the non-resonant hy- 
perfine state. However, at certain magnetic fields, the 
hyperfine level lifetime plummets, leading to a decrease 
in the hole amplitude and the complete extinction of any 
antiholes. As we will show, these points of holeburning 
extinction occur when the magnetic field brings the Tm 
hyperfine transition into resonance with an Al transition. 

Extinction points can be seen in Figure[lj which shows 
the holeburning structure of the 3 H 6 (0) —UH 4 (0) transi¬ 
tion of Tm 3+ :YAG as a function of magnetic field magni¬ 
tude. The experimental procedure used to generate the 
spectrum in Figure |T] was as follows. A spectral hole was 
burnt with a weak, 100 ms pulse, and after a delay of 
10 ms, a chirped pulse was used to read out the resulting 
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FIG. 1: (Color online) Holeburning transmission spectrum 
in Tm 3+ :YAG as a function of magnetic field along [001]. The 
color scale is logarithmic. One sidehole, and two antiholes are 
visible on each side of the central hole. These are labelled by 
their splitting from the central hole. Near complete extinc¬ 
tion of the spectral hole occurs at 3.1 mT and 6.2 mT. Each 
horizontal spectrum is an average of eight shots. Between 
each shot the magnetic field was switched off, or to the near¬ 
est extinction point, to prevent accumulation of the spectral 
hole. 


holeburning structure. Each horizontal slice in the fig¬ 
ure is the holeburning spectrum for a particular value of 
the applied magnetic field. The magnetic field was along 
the [001] direction of the crystal. There are six magnet¬ 
ically inequivalent Tm sites in Tm 3+ :YAG, labelled as 
sites 1 to 6 (see Ref. USD, and for a field orientation of 
[001] four sites (sites 3 to 6) are equivalent. To eliminate 
the other two sites, sites 1 and 2, from the hole burn¬ 
ing spectrum, the laser beam direction was chosen along 
[110] and the polarization along [001]. The holeburning 
spectrum, therefore, shows the structure of a single class 
of ions: one sidehole at A e , the excited state splitting, as 
well as two antiholes at A g , the ground state splitting, 
and A g — A e . A third antihole is expected at A g + A e , 
but for the ma gne tic field orientation chosen, this is too 
weak to be seerpSl. 

The figure shows that when the magnetic field is in¬ 
creased from zero the central hole amplitude increases 
as the ground state hyperfine levels split and the hy¬ 
perfine lifetime starts to rise. However, at two fields, 
3.1 mT and 6.2 mT, the central hole amplitude suddenly 
drops, and the antiholes disappear, indicating that the 
hyperfine ground state lifetime has decreased to less than 
the metastable state lifetime. The lifetime is substan¬ 
tially shortened over a wide region, 0.6 mT, which corre- 
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FIG. 2: (Color online) Holeburning extinction due to res¬ 
onances between the Tm ground state splitting and the A1 
quadrupole splitting. The A1 level structure is shown in the 
inset, and the horizontal red lines show the splittings between 
these levels. The dotted red horizontal line indicates an A1 
transition (±1/2 —> ±5/2) that is extremely weak. Each time 
one of the A1 splittings is in resonance with the Tm splitting 
(blue line), an extinction point can be expected. In Figure [I] 
we observe the resonances at B ex t and 2 B ex t. 


sponds to a range in liyperfine splitting of approximately 
200 kHz. 

The hole extinction behavior can be explained by a 
resonant cross relaxation between Tm and Al. The in¬ 
teraction leads to a mutual spin flip-flop whereby the Tm 
excitation is transferred to the Al, filling in the Tm spec¬ 
tral hole. In Figure [T] we see two extinction points, with 
the second at twice the held of the first. This is clear evi¬ 
dence that the interacting species is Al. Al has a nuclear 
spin of / = 5/2 and therefore a nuclear quadrupole mo¬ 
ment. In the high symmetry Al sites of YAG, this tensor 
is cylindrically symmetric, giving rise to three doubly de¬ 
generate hyperhne levels in zero Held with splittings A q 
and 2 Aq. Al also has an isotropic nuclear Zeeman ten¬ 
sor, which lifts the double degeneracy of the three sets of 
levels in a magnetic Held. However, because the nuclear 
magnetic moment of Al is more than an order of mag¬ 
nitude smaller than the Tm enhanced nuclear moment 
it has little effect on the position or behavior of the ex¬ 
tinction points, and we will ignore it for now. It will be 
considered in detail in Section HV Bl 

With this assumption, the energy level splittings of the 
Tm and Al ions can be drawn as in Figure [2] The two 
Tm levels split linearly with magnetic field, while the Al 
levels are constant. Extinction points can occur when 
the Tm ±1/2 —> —1/2 transition is resonant with one of 
the Al transitions. The two extinction points in Figure [l] 
correspond to the first two Al transitions, the ±1/2 —> 


±3/2 transition at A q and the ±3/2 —> ±5/2 transition 
at 2 Aq. A third extinction point is expected at 3 Aq 
when the Tm splitting crosses the third Al transition 
±1/2 —> ±5/2, although this transition is very weak so 
this point is likely to have less effect on the Tm lifetime 
than the other two extinction points. 

The Al quadrupole splitting Aq can be written 


A Q = 


6 Q p 

4/(27- 1) 


(1) 


where Q p is the quadrupole parameter. Al occupies 
two sites in YAG with very different quadrupole split¬ 
tings: an octahedral (C^ symmetry) site with Q p = 
0.632 MHz and a tetrahedral site (S 4 symmetry) with 
Q p = 6.017 MHi® J . The extinction points are caused 
by cross relaxation with the tetrahedral site, which has 
A q = 0.90 MHz. This is in good agreement with the 
experimental position of the extinction point at a Tm 
splitting of 0.91 ± 0.02 MHz, measured from the hole¬ 
burning spectrum in Figure [l] 

Extinction points like those in Figure [l] can be seen for 
any direction of the magnetic field. They occur at fixed 
Tm ground state splittings, equal to the Al splittings. 
However, because the nuclear Zeeman tensor of Tm in 
YAG is highly anisotropic, the magnetic field values at 
which the extinction points occur vary with the field di¬ 
rection. Since one component of the Tm ground state 
Zeeman tensor (q.y « 403 MHz/T) is more than twenty 
times larger than the other two components, the ground 
state splitting can be approximated from the projection 
of this y component onto the magnetic field direction. For 
the [001] orientation of the magnetic field in Figure]!] the 
four sites contributing to the structure each have their y 
axis 45° from the magnetic field direction. The Zeeman 
splitting rate of the ground state is therefore = ^,or 
285 MHz/T. This gives splittings of 0.88 and 1.77 MHz 
at the two extinction points, in agreement with the ob¬ 
served splittings. 

Figure [3] shows the holeburning structure for a differ¬ 
ent magnetic field direction, [110]. In this direction, as 
for [001], sites 3 to 6 are equivalent. The laser beam ori¬ 
entation was again chosen to eliminate site 2 but for this 
case, the polarization was chosen along [ 110 ], which al¬ 
lows site 1 to be seen. The spectrum, therefore, displays 
two sets of sideholes and antiholes. The inner antiholes 
belong to sites 3 to 6 while the outer antiholes belong 
to site 1. For this magnetic field direction the strongest 
antihole, at A g — A e , is visible for both sites 3 to 6 and 
site 1 while one of the weaker antiholes, at A g , is only 
visible for sites 3 to 6 . 

The two sets of antiholes display different extinction 
behavior. The outer, site 1 antiholes disappear at 2.25 
and 4.5 mT, while the inner, sites 3 to 6 , antiholes disap¬ 
pear at 4.5 and 8.8 mT. Therefore, the extinction point 
at 4.5 mT is an extinction point for both types of sites, 
while the 2.25 mT point is only an extinction point for 
site 1 and the 8.8 mT one only for sites 3 to 6 . This be¬ 
havior can be understood by considering the splitting of 
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FIG. 3: (Color online) Holeburning transmission spectrum in 
Tm 3+ :YAG as a function of magnetic field along [110]. The 
color scale is logarithmic. The inset is a magnified version of 
the region inside the rectangle: the antihole structure around 
the first extinction point. A clear inversion of the antihole 
is seen, which is explained in Section [IV A| Horizontal slices 
below 7.1 mT are 8-shot averages as described in the caption 
of Figure[Tj Extinction of the spectral hole occurs at 2.25 mT 
and 4.5 mT, 6.75 mT and sa 8.8 mT. Above 7.1 mT, the mag¬ 
netic field could not be automatically switched off between 
shots to erase the structure, so each spectrum is a single shot, 
and the magnetic field was switched off manually between 
each spectrum. The additional features seen in these spectra 
are artifacts of this manual method. 

these two sets of sites. The y axis of the Zeeman tensor 
for site 1 is directed along [110], the magnetic field direc¬ 
tion, so the ground state splitting rate takes its maximal 
value = 7 , y . Meanwhile, sites 3 to 6 have their y 
direction 60° from the magnetic field, so for these sites 
= Therefore, at 2.25 mT, site 1 is resonant with 
the lower A1 splitting Aq, while at 4.5 mT, site 1 is now 
resonant with the higher splitting 2 Aq and sites 3 to 6, 
which have a ground state splitting exactly half that of 
site 1, are resonant with the lower splitting A q. Finally, 
at S3 8.8 mT sites 3 to 6 are resonant with 2 Aq. 

Careful examination of Figure [3] also shows that the 
amplitude of the site 1 antihole decreases at 6.75 mT, to 
approximately 30 % of its original size. At this magnetic 
field, the Tm ground state splitting is 2.7 MHz, resonant 
with the weak 3 Aq, ± 3 —> ±| transition. 

Two things should be noted from the holeburning spec¬ 
tra of Figures [l] and [3] First, there is no visible deviation 
of the position of the antiholes near the resonance points. 
This indicates that the interaction is sufficiently weak 
that it does not substantially perturb the energy level 
structure of the Tm and A1 ions near resonance. Second, 


the holeburning extinction occurs over a very large range 
of Tm transition frequencies, of the order of 200 kHz. 
Since the interaction itself is too weak to cause such a 
wide extinction range, this demonstrates that other fac¬ 
tors, in particular the Tm inhomogeneous broadening of 
0(100) kHz, influence the holeburning structure near the 
resonance. 


III. MEASUREMENT OF HYPERFINE 
LIFETIMES USING STIMULATED ECHOES 

In the previous section, we showed that the hyperfine 
lifetime is reduced at points of resonance between Tm 
and Al, resulting in a reduction in the central hole and 
antihole amplitude. The value of the lifetime, at and 
away from resonance, could be determined from a hole¬ 
burning spectrum by measuring the hole amplitude as a 
function of the delay between the holeburning pulse and 
the readout chirp. However, in a spectral holeburning 
sequence, the distribution of ions amongst the states of 
the system at the end of the holeburning pulse is highly 
dependent on the hyperfine lifetime itself, meaning that 
ions at resonance and those away from resonance are pre¬ 
pared in completely different ways by the holeburning. 
This makes analyzing the hole decay somewhat compli¬ 
cated at this stage. An easier method is to use a stimu¬ 
lated (3-pulse) echo sequence, and analyze the echo am¬ 
plitude as a function of the delay between the second and 
third pulses. In this sequence the first two pulses prepare 
the system in a consistent state, greatly simplifying the 
extraction of the hyperfine lifetime. In this section, we 
present the theory and experimental data for stimulated 
echoes at, and away from, the resonance. 


A. Relaxation dynamics 

To determine the hyperfine lifetime from a stimulated 
echo decay in a multilevel system like Tm 3+ :YAG, it is 
necessary to have a model for the shape of the decay 
curve as a function of the lifetimes and branching ratios 
of each of the levels involved in the decay process. For 
Tm 3 +:YAG, a rate equations model with four levels - 
two hyperfine ground states, one excited state, and a 
metastable state - is sufficient. We are interested in the 
decay behavior so we consider a set of rate equations with 
no driving terms. This four level system is illustrated in 
Figure [4] We will start by assuming that the decay rates 
are the same for all atoms in the system, which applies 
when the atoms are far from an extinction point. We 
will consider the situation when this assumption does 
not apply later. 
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FIG. 4: (Color online) 4-level system used in the rate equa¬ 
tion model. R 0 , R m and R g are the decay rates of the optical, 
metastable, and ground states respectively while p 0 , p m and 
p g are the branching ratios that describe how each of these 
levels decays. 


The rate equations for the system are: 
cIti R, 

= -y («2 - n i) +p 0 ( 1 -p g )R 0 n 3 + (1 -Pm)Rm,ni 
dri2 R q , , _ _ 

-jjT = - n i) +PoPgR 0 n 3 +PmR m ni 


7, — (1 Po)R 0 n 3 Rm 7/4 

at 

( 2 ) 

Where the R parameters are transition rates and the p 
parameters branching ratios, as defined in Figure[4j The 
solution to this linear system of equations is a sum of 
four exponentials with decay rates of 0, R g , R m and R 0 - 
The explicit solution is given in Appendix [A] 

We now apply these rate equations to a stimulated 
echo experiment. The stimulated echo pulse sequence 
used in the experiment involves three co-linear pulses, 
with the first and second pulses separated by time r short 
compared to the coherence time X 2 , and the second and 
third pulses separated by a time T w much longer than 
X 2 . We consider that the pulses are resonant with the 
transition |2) —I 13), and the initial population is evenly 
distributed between the two ground states, rq = 712 = 
0.5. The first two pulses engrave a spectral population 
grating between the levels |2) and |3) with a period given 
by the pulse spacing r. The populations of the system 
can be described by a four element population vector 
n(t, v). After the engraving step, 

n(T,v)=n eq + £(T,v)ri gr (3) 

where n eq = [f, 5 , 0 , 0 ] describes the equilibrium popu¬ 
lation distribution, n gr = [ 0 , - 1 , 1 , 0 ] describes the popu¬ 
lation inversion caused by the engraving pulses, while the 
frequency dependence of the grating is given by £(t,v). 


In the simplest case, this is of the form 

£.{r,v) = ^0 (1 + cos (27n/r + 0)) (4) 

for frequencies within the grating bandwidth. £ 0 <C 1 is 
the grating depth. 

During the wait time between the second and third 
pulses, the grating will degrade. The evolution of n(t, u) 
after the engraving step can be written as 

n{t,v) = U(t — r)n(T,v) (5) 

where the evolution matrix U can be derived from the 
solution to the rate equations (Equations <§)• Since the 
evolution matrix does not affect the equilibrium popula¬ 
tion n eq , the population can then be written 

n(t, v) = n eq + £(t, v)U [t — T)n gr (6) 

Thus, the evolution of the grating is described simply by 
the evolution of n gr under U. 

The third pulse of the stimulated echo sequence scat¬ 
ters off the spectral grating, forming an echo at time r 
after the third pulse. The echo amplitude is proportional 
to the visibility of the spectral grating at the time of the 
third pulse, and can be written 

A(T W ) cx n T gr U{T w )n gr (7) 

The echo intensity, which is the parameter measured in 
the experiment, is then 

I(T W ) =C 0 (nl r U(T w )n gr ) 2 (8) 

where Cq is a proportionality constant. Since U is de¬ 
rived from the rate equations, there are seven unknowns 
in this equation: Co, the three decay rates R 0 , R m , and 
R g , and the three branching ratios p Q , p m and p g . Of 
these, p g , which is the branching ratio from the excited 
state to the initial hyperfine ground state, can be fixed at 
p g = 1. In general, p g is the same as the branching ratio 
of excitation from the hyperfine ground states to the ex¬ 
cited state if the decay is direct or if the time spent in any 
intermediate levels is short, which is the situation that 
applies to Tm 3+ :YAG. Due to the highly anisotropic Zee- 
man interaction in Tm 3+ :YAG, the excitation branching 
ratio is close to one along all but a few special directions. 

B. Decay away from resonance 

Equation (JsJ) gives a model of the echo decay curve as a 
function of six parameters. To extract these parameters 
we measured the stimulated echo decay at a magnetic 
field away from any extinction points. All measurements 
were performed in a helium bath cryostat with the sam¬ 
ple temperature maintained below 2.6 K. Close control of 
the temperature is necessary for Tm 3+ :YAG as the hy¬ 
perfine ground state lifetime is temperature dependent 
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FIG. 5: (Color online) Stimulated echo intensity decay for 
magnetic fields along the [001] direction far from resonance 
(4.56 mT, blue crosses), and close to resonance (6.2 mT, green 
dots). T w is the delay between the second and third pulses of 
the sequence. The red line is a fit to the far from resonance 
data using Equation (|8|, while the purple line is a fit to the 
near resonance data using Equation HD- 


above approximately 2.8 K. An external coil was used to 
supply a magnetic field of 4.56 mT parallel to the [001] 
direction. The laser beam was oriented along [110] and 
polarized along [001], which, as noted in the previous 
section, means that only the four equivalent sites 3 to 6 
interact with the laser. 

The laser power was 1.6 mW, focused onto the crystal. 
In the first part of the stimulated echo sequence, two 

1 ps square pulses separated by 5 ps were used to create 
the population grating, which was read out with a third, 
identical pulse at delays between 100 ps and 60 ms. The 
repetition rate was slow, 0.3 Hz, and after every shot the 
magnetic field was switched to the extinction point at 
3.1 mT for 200 ms to erase the population grating. 

The decay of the stimulated echo intensity is shown in 
Figure [5] Also shown is a fit to the data using the model 
of Equation (§. Three decay regions can be identified in 
the figure. At first, T w < 2 ms, the curve is dominated 
by the rapid excited state decay, with a rate R 0 = (2.05± 
0.04) x 10 3 s _1 (T° = 0.49 ± 0.01 ms). For intermediate 

2 < T w < 30 ms, the metastable state decay, with R m = 
94 ± 4 s _1 (T™ = 10.6 ± 0.4 ms) dominates, while at 
larger times, the slow ground state decay, with R g = 4.4± 
0.8 s~ 3 (Tf = 230 ±40 ms) begins to become important. 
The fit to the data also gives the branching ratios p Q = 
0.30 ± 0.02 and p m = 0.67 ± 0.01. As stated above, p g is 


fixed at 1. 

The fitted values for the lifetime of the optical 
and metastable levels, 0.49 ms and 10.6 ms respec¬ 
tively, a re in good agreement with previously published 
values 12 20 21 .The hyperfine level lifetime is considerably 
shorter than the lifetimes of several seconds that can be 
seen in holeburning measurements at similar fields. This 
suggests that some degree of spectral diffusion is con¬ 
tributing to the decay process. We chose not to add 
spectral diffusion to our model, because the only param¬ 
eter affected by this inclusion is the hyperfine lifetime, 
and it was not a priority of our experiment to measure an 
accurate hyperfine lifetime in the non-resonant regime. 


C. Decay near resonance 

The fitted parameters gained in the previous section 
for the decay away from resonance can be used to build 
a model for the decay at resonance. The crucial difference 
between these two situations is that at resonance, the hy¬ 
perfine state decay can no longer be described by a single 
decay rate. This occurs because the Tm hyperfine tran¬ 
sition is inhomogeneously broadened, so that the lifetime 
of each ion depends on how far away it is from resonance 
with the A1 transition. The ensemble must be described 
by a distribution of hyperfine lifetimes. The A1 transi¬ 
tion is also inhomogeneously broadened, but because the 
A1 magnetic moment is small, this broadening is much 
smaller than the Tm broadening and will be ignored. 

To model the effect of the Tm inhomogeneous broad¬ 
ening on the echo near resonance, we first write the hy¬ 
perfine decay rate R g as a function of the distance from 
resonance with the A1 ion: 

R g (v) = r 0 + (r res - r 0 )h{v - v A i) (9) 

where tq is the decay rate far from resonance measured in 
the previous section, r res is the rate exactly on resonance, 
and hiy—VAi) is the lineshape of the resonant interaction 
between A1 and Tm, centered on the A1 splitting vai- 
For the moment, we will assume this function is a single 
Lorentzian with an amplitude of one. In Section [I V B | we 
will examine h. in more detail. The width of h is T/y, 
the flip-flop resonance width, which is dependent on the 
strength of the Al-Tm interaction. 

As the hyperfine level decay rate is now a function of 
the Tm ground state splitting frequency, the contribu¬ 
tion to the echo amplitude for the subset of ions with 
hyperfine frequency v is 

A(T, u, A g ) = Cg /2 ( nl r U(T)n gr ) - A g ) (10) 

where G{v — A g ) is the inhomogeneously broadened line- 
shape of the hyperfine transition, centered on the fre¬ 
quency A g (B). The echo intensity, then, is the square of 
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the integral of A(T, v, A g ) over the lineshape G{v - A g ): IV. MODELLING THE HOLEBURNING 

STRUCTURE 


I(T,A g )=C 0 



(■ n J gr U(T)n gr ) G{v 



The maximum effect on the echo is seen when the Trn 
and A1 transitions are in resonance, A g = vai- At 
this point, the overlap between the flip-flop resonance 
lineshape h(v — vai) and the inhomogeneous lineshape 
G{y — A g ) is greatest. 


Equation (111 is dependent on the parameters deter¬ 
mined in the previous section, as well as three new pa¬ 
rameters: the maximum hyperfine decay rate r res , the 
flip-flop resonance width Tf f and the hyperfine inhomo¬ 
geneous broadening, T^. This last parameter can be 
determined from the holeburning spectrum in Figure [T] 
The lineshape of the antiholes in this spectrum is a con¬ 
volution of the central hole lineshape with the hyperfine 
inhomogeneous broadening function G. A fit to the anti¬ 
hole lineshape gives G as a Gaussian function with T^ 
of 120 kHz (FWHM). 


In the previous section, we showed that the stimulated 
echo decay near resonance could be well explained by a 
model that took into account the inhomogeneous broad¬ 
ening of the Tm hyperfine transition. One of the implicit 
assumptions in this model was that the inhomogeneous 
broadening is much greater than the Tm-Al flip flop res¬ 
onance width. In the present section, we substantiate 
this assumption by investigating the holeburning struc¬ 
ture near resonance. 

As previously mentioned, the holeburning structure is 
slightly difficult to model because both the preparation 
and readout steps of the experiment are highly depen¬ 
dent on the hyperfine lifetime, unlike the echo decay, for 
which only the readout step is lifetime dependent. To ac¬ 
curately model the holeburning structure, therefore, we 
make use of the lifetimes and the branching ratios cal¬ 
culated in the previous section from the stimulated echo 
decays. We will concentrate on the anti-hole structure, 
as, for reasons we will shortly explain, this is more sensi¬ 
tive to the details of the resonant Tm-Al interaction than 
either the hole structure or the echo decay. 


The other two parameters, r res and T^y, cannot be 
determined solely from a fit to the echo decay curve at 
resonance because they both have a similar effect on the 
curve. However, the combination of the echo decay fit 
and the fit to the antihole structure near resonance (Sec¬ 
tion IV C) determine that r res = (2 ± 0.8) x 10 3 s _1 and 
T ff = 12 ± 3 kHz. 


Figure [5] shows the echo decay for the 1.8 MHz extinc¬ 
tion point at 6.2 mT along the [001] direction. The echo 
decays much more rapidly than the far-from-resonance 
curve shown in the same figure, and has only two clear de¬ 
cay regions: at short times, the decay due to the excited 
state and at longer times, that due to the metastable 
state. The purple line in the graph shows the model pre¬ 
diction of Equation mi using T inh = 120 kHz, r res = 
2 x 10 3 s _1 (T-f(res) = 0.5 ms), Tff = 12 kHz, and 
the decay rates and branching ratios from the far-from- 
resonance fit. This calculated Tf for the hyperfine level is 
substantially shorter than the metastable state lifetime, 
and of the order of the excited state lifetime. 

In Figure [5] the model of Equation mi matches the 
experimental decay curve well. There is a slight discrep¬ 
ancy at medium to long times, which can be attributed to 
the extreme simplicity of the model used. In particular, 
it was assumed that the effect on the hyperfine lifetime 
of the resonance with the A1 could be represented by a 
single Lorentzian fixed at the center of the resonance. 
This assumption only holds if the A1 Zeeman structure 
is ignored, and we show in Section [TV B| that accurately 
modelling the behavior near resonance requires taking 
the Zeeman structure into account. 


A. Antihole and hole structure 

To model the hole and antihole structure, we use the 
four level system of Figure [4] where the driving field is 
again tuned to the |2) -A |3) transition. We assume that 
the holeburning step is not saturated, so that the shape 
of the spectral hole C ( v — U>) centered around the opti¬ 
cal frequency v 0 is invariant with the distance from reso¬ 
nance A g — vai- Then, only the amplitude of the spectral 
hole changes with A g — vai ■ We begin with the structure 
of the central hole and the anti-hole for ions with a fixed 
ground state splitting A g . 

The central hole amplitude for a particular ground 
state splitting is dependent on the population 772 at the 
end of the holeburning and wait steps. The population 
after the holeburning step 712 ( 4 , A g — vai) can be cal¬ 
culated from a rate equations model similar to Equa¬ 
tions ([ 2 ]) with a driving term added, in which the ground 
state relaxation rate is allowed to vary as in Equation (|9j) . 
Then, the population after the wait step 712 (t w , A g — vai) 
can be determined by applying the solution to Equations 
([ 2 ]) with the same varying R g to 772 ( 4 , A g — vai)- The 
logarithm of the hole amplitude is then proportional to 

^/i(Ag vai) = weq ^12(4,, A g vai) ( 12 ) 

where n eq = 0.5 is the equilibrium population of the level. 
A similar expression exists for the antihole: 

k ah {A g - vai) = n eq - ni(t W: A g - v A i) (13) 

Equation [12] gives the contribution of ions at a par¬ 
ticular frequency A g to the hole amplitude. If we now 






take into account the inhomogeneous broadening of the 
hyperfine ground state, the hole structure for an average 
ground state splitting of A g is given by integrating kh 
over the inhomogeneous lineshape 


H h (u, A g ) =C,{v- z/ 0 ) 


/• + OO 


k h {Ag + a - VAi)G(a)da 

J — OO 

(14) 

The integral is simply a convolution of kh and G. This 
convolution means that the hole amplitude is reduced 
over a range of Trn hyperfine splittings comparable to the 
inhomogeneous broadening, a range much bigger than the 
width of kh itself. This explains why the experimental 
holeburning spectra in Figures [l] and [3] show holeburning 
extinction over such large regions. A second important 
consequence of the convolution is that any structure in 
k h smaller than the inhomogeneous broadening will not 
be visible in Hh- 

In contrast to the central hole, the anti-holes are sensi¬ 
tive to the structure of the resonance. To prove this, we 
consider the antihole at A g as the equations are some¬ 
what simpler for this case. We present results for the 
antihole at A g — A e , the one most commonly seen in the 
experiment, in the Appendix [B] 

The antihole centered at an optical frequency of A g 
can be described by 

+OO 

H ah (v,Ag) = - J k a h{5v+a-VAi)G(5v+a-A g )C,(a)da 

— OO 

(15) 

where 8v = v — v o is the optical detuning. In this case, 
the convolution is between the hole lineshape and the 
product of k a h(8 v — vai) and G(8v — A g ). Because the 
functions k a h and G are multiplied and not convoluted as 
for the structure of the hole (Equation (fl4|), the antihole 
near the resonance will be sensitive to any structure in 
k a h larger than the width of £. 

Equation (15) predicts that in a two-dimensional hole¬ 
burning spectrum such as Figure [3] the resonance point 
will appear as a vertical line of reduced antihole ampli¬ 
tude which is separated from the central hole by the A1 
splitting and whose vertical extent is given by the in¬ 
homogeneous broadening. This broadly agrees with the 
appearance of antiholes in the figure. 

Many of the extinction points in Figure [3] show not 
just reduction but inversion of the antihole. This is most 
clearly seen in the inset, a magnified view of the first 
extinction point at 2.25 mT. The inversion can be ex¬ 
plained by considering the dynamics of the holeburning 
process. At resonance, we have shown that the hyperfine 
state lifetime is shorter than the metastable state life¬ 
time. Importantly, it is also shorter than the holeburn¬ 
ing pulse length, meaning that the two hyperfine levels 
can come into equilibrium during the holeburning pulse. 
As the resonant level is continually being emptied by the 
laser, the relaxation between the two hyperfine levels will 
begin to empty the non-resonant one. Therefore, at res¬ 
onance the holeburning process pumps population out 


of both hyperfine levels to be stored in the metastable 
state, burning a sidehole and not an antihole in the non¬ 
resonant level. 

Another aspect of the antihole structure that is visible 
in the inset of Figure [3] is that the inverted region is diag¬ 
onal, and not vertical. This diagonal shape indicates that 
the antihole extinction function k a h is structured, made 
up of multiple resonances whose separations are larger 
than the individual resonance width. The diagonal in¬ 
verted structure, therefore, is made up of a number of off¬ 
set vertical resonances as described by the Equation [15] 
Multiple closely spaced resonances arise in Tm 3+ :YAG 
when we take into account the Zeeman structure of Al, 
which is the topic of the next section. 


B. Role of the Al Zeeman splitting 

Thus far we have ignored the effect of the small Al 
Zeeman splitting in order to simplify the models for the 
holeburning structure and stimulated echo decay. We saw 
in the previous section that the Al Zeeman splitting does 
have an effect on the antihole structure near resonance, 
so in this section we include it into the model for the 
antihole structure. 

Ignoring the Al Zeeman interaction substantially sim¬ 
plified the models of the previous sections because it 
is equivalent to assuming that all Al ions in the crys¬ 
tal have the same hyperfine structure. When we take 
the Al Zeeman tensor into account, this is no longer 
the case. While the Al Zeeman tensor, which has 
7 ai = 11.1 MHz/T, is isotropic, when combined with 
the anisotropic quadrupole tensor it leads to energy lev¬ 
els splittings that are dependent on the orientation of 
the quadrupole axis with respect to the applied magnetic 
field. This means that Al ions with different orientations 
will have different energy level structure, and therefore 
lead to differently structured extinction points when they 
interact with Tm. 

The tetrahedral Al ions in YAG have three different 
orientations for the quadrupole axis, along the three 
(001) directionj2H. For magnetic fields along [001] and 
[110], there are, therefore, three possible orientations of 
the axis relative to the field: parallel, perpendicular, and 
at an angle of 45°. Each will lead to a different pattern 
of splitting of the Al hyperfine levels. As an example, 
we show the splitting for the case of 45° in Figure [6] 
The magnetic field lifts the double degeneracy of all three 
quadrupole levels. A similar pattern with larger splitting 
is seen for a field parallel to the quadrupole axis, while 
for a field perpendicular to the axis only the m z = 
levels split. 

Since the Al splittings are now dependent on the ori¬ 
entation of the Al sites, the extinction point structure 
will be dependent on which Al ions a particular Tm site 
interacts with. This can be determined from the crystal 
structure. Tm in YAG occupies six magnetically inequiv¬ 
alent sites of D 2 symmetry (3 mutually perpendicular C 2 
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Magnetic field (mT) 

FIG. 6: Hyperfine transition frequencies of A1 in YAG for a 
magnetic field 45° from the direction of the quadrupole quan¬ 
tization axis. Solid lines indicate strong transitions (transi¬ 
tion probability above 0.1) while dashed lines indicate weak 
transitions (probability below 0.01). 




Splitting frequency (MHz) 


symmetry axes). The difference between these six sites 
is the orientation of their symmetry axes. Each pair of 
sites (1 and 2, 3 and 4, 5 and 6) has one C 2 axis along 
a (001) direction, and the other two along perpendicular 
(110) directions. We will denote the (001) direction of a 
site as the polar axis. 

Each Tm ion has two nearest neighbor tetrahedral A1 
ions at 3.0 A in either direction along the polar axi d 23 * 24 ( 
The quadrupole axis for these two A1 ions is also along 
the polar axis. Thus, for site 1 and 2 Tm ions the nearest 
neighbor A1 ions have their quadrupole axis along [001], 
for sites 3 and 4 along [100] and for sites 5 and 6 along 
[ 010 ], 

The next nearest neighbor tetrahedral A1 ions are 4 
ions at a distance of 3.6 A. These 4 ions divide into two 
pairs with quadrupole axes along the two (001) directions 
perpendicular to the polar axis. The third nearest neigh¬ 
bors are considerably further away than the next nearest 
neighbors, at 5.6 A, and will not be considered here. 

To show how this information can be used to determine 
the A1 resonance frequencies for a particular Tm site, we 
shall use the extinction point at 2.25 mT in Figure [3] as 
an example. In this figure, the magnetic field is along 
[110] and the extinction point is due to ions in site 1, 
whose polar axis is [001]. Thus, their nearest neighbour 
A1 ions will have their quadrupole axis along [001] as 
well, perpendicular to the field, while the next nearest 
neighbours have their quadrupole axis at 45° to the field. 


FIG. 7: (Color online) (a) Site 1 antihole amplitude as a 
function of the Tm ground state splitting for a magnetic field 
along [110]. All three A1 resonances, at 0.9, 1.8 and 2.7 MHz, 
are visible, (b) Position of the resonant A1 transitions for the 
two nearest neighbor ions to site 1, whose polar axis is per¬ 
pendicular to the field. The probabilities for the transitions 
above 2.5 MHz are multiplied by 100 to make them visible, (c) 
Position of resonant A1 transitions for the four next nearest 
neighbor ions ( quadrupole axis 45° to the field). 


Figure [7] compares the energy level splitting for these 
sets of A1 ions (bottom two plots) to the depth of the 
antihole as a function of magnetic field, top, which was 
obtained by cutting diagonally across Figure [3] through 
the center of the antihole. In this figure, vertical lines in¬ 
dicate calculated A1 hyperfine transition frequencies for 
the nearest and next nearest neighbour sites. Their am¬ 
plitudes are given by the corresponding transition proba¬ 
bility, which can be calculated for an RF transition from 
the expectation value of the magnetic dipole operator. 

The positions and widths of the extinction points in 
Figure[?]correlate well with the A1 structure: the first two 
extinction points are relatively narrow due to the small 
A1 splittings, while the third, weaker extinction point is 
much broader due to the large spread in A1 structure. 
The degree of inversion of the antihole, which decreases 
with increasing splitting, also agrees with the general de- 
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cline in A1 transition probabilities with increasing split¬ 
ting. 


C. Comparison between model and antihole 
structure 


The calculation of the A1 Zeeman structure in the pre¬ 
vious section can be used to more accurately model the 
Tm antihole structure near resonance. To modify the 
antihole structure model of Section IIVI to include the A1 
structure, it is sufficient to feed an altered function for 
R g into the rate equation models used: 


R g (v) =r 0 + (r res - r 0 ) 




(16) 


where vai,% are the different A1 transition frequencies that 
contribute to the resonance. 

We will again apply this to the site 1 extinction point 
for a 2.25 mT magnetic field along [110]. This is the 
simplest case to model as only one site contributes to the 
antihole. Since the visible antihole is at A g — A e , we use 
the antihole structure equation suitable for this antihole 


(Equation Bl). We choose to only consider the effect 
of nearest neighbour A1 sites as these should dominate 
the interaction. The nearest neighbor ions lead to two 
resonances vai at 0.86 and 0.94 MHz (see Figure]?]). The 
lineshape of the hyperfine decay rate R g (v) is then given 
by Equation (161. 


The model is dependent on the inhomogeneous broad¬ 
ening lineshape G(5v— ( A g — A e )) and the hole lineshape 
£(z/), which can be measured from the holeburning spec¬ 
trum away from resonance. As previously mentioned, G 
is a Gaussian function with FWHM 120 kHz, while £ is 
a Lorentzian with a lincwidth of 50 kHz. 

The antihole structure is also dependent on r res and 
Tff, parameters that occur in the model for the stimu¬ 
lated echo decay shape of Section |mc| As described in 
that section, the values of r res and Tf f that best fit both 
the echo decay and the antihole structure near resonance 
are r res = (2 ± 0.8) x 10 3 s _1 and Tff = 12 ± 3 kHz. 

Using the parameter values described above, Figure [8] 
compares the model spectrum to the experimental data 
for several different values of the ground state splitting 
near the resonance. The model matches the general be¬ 
havior well, showing first a partial inversion of the high 
frequency components of the antihole as these ions come 
into resonance (0.82 MHz), then a centralized, complete 
inversion (0.91 MHz), and finally a partial inversion on 
the low frequency side (0.99 MHz). The amplitude of the 
antihole far from resonance (0.57 MHz), and the ampli¬ 
tude of the inverted peak are well matched by the model, 
although the experimental spectrum takes longer to re¬ 
turn to the original antihole size, suggesting that the 
wings of the resonance are slightly wider than the model 
predicts. This could be due to a contribution from the 
next-nearest neighbors, which were ignored in the model. 



FIG. 8: (Color online) Comparison between experimental an¬ 
tihole structure (blue solid lines) and modeled structure (red 
dashed lines) near the site 1 resonance at 2.25 mT along [110]. 
The different traces are different Tm ground state splittings, 
with the value of the splitting listed for each trace in MHz. 
The resonance is at approximately 0.91 MHz. 



Tm ground state splitting (MHz) 


FIG. 9: (Color online) Central hole amplitude as a function 
of the Tm ground state splitting for the 2.25 mT extinction 
point of site 1 along [110]. The experimental spectrum is 
given by the blue dots while the model of Equation ( |14[ ) is 
shown by the red line. 
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The same set of parameters as used above can be used 
to model the central hole structure near the resonance, 
based on Equation (141. Figure [9] shows the experimen¬ 
tal amplitude of the central hole compared to that of the 
model as a function of the Trn ground state splitting. 
The resonance is relatively shallow because half the signal 
comes from ions in sites 3 to 6, which are not at an extinc¬ 
tion point. The model is a close fit to the experimental 
data in the region of the resonance. Outside the reso¬ 
nance region, the experimental amplitude increases with 
the splitting, an indication that the lifetime is length¬ 
ening. This shows that there are other factors affecting 
the lifetime apart from the Al-Tm cross-relaxation reso¬ 
nances. 


V. DISCUSSION 

We have shown that, when the Tm and A1 nuclear 
splittings coincide, relaxation of the Tm spin happens 
rapidly through a resonant flip-flop interaction with 
nearby A1 ions. The relaxation is completed within a few 
milliseconds, meaning that A1 thermalization rate must 
be consistent with that time scale. 

The Tm-Al flip-flop is expected to be driven by the nu¬ 
clear magnetic dipole-dipole interaction. This interaction 
can be modeled using the well-known crystal structure 
of YAG and the Tm and A1 gyromagnetic tensors. Un¬ 
derstanding the relaxation process requires, in addition, 
some model for the stochastic flipping of the A1 spins, 
which results from their coupling to each other. The in¬ 
teractions between A1 ions are likely to be affected by 
the existence of a small frozen core around the Tm ion 
arising from the enhanced nuclear magnetic moment of 
Tm. This frozen core complicates the model of the A1 
spin bath dynamics, as it means that A1 near the Tm 
relax at different rates to those further away. We defer 
the theoretical analysis of these spin bath dynamics, and 
the resulting Tm-Al flip-flop rate, to a future paper. 

For the practical applications of Tm 3+ :YAG to quan¬ 
tum and classical information processing, the reso¬ 
nant Tm-Al interactions studied here could prove use¬ 
ful for manipulating the Tm hyperfine level populations. 
Switching the magnetic field to an extinction point pro¬ 
vides a way of equalising the populations of the two levels 
without reducing the field completely to zero. It is also 
possible to manipulate the hyperfine states of the dif¬ 
ferent Tm sites independently, as they have extinction 
points at different fields. 


VI. CONCLUSION 

At well defined magnetic fields < 10 mT, the lifetime 
of the Tm 3+ :YAG hyperfine ground states decreases sub¬ 
stantially, which is seen as a reduction in the holeburn¬ 
ing amplitude. This was attributed to a resonant cross¬ 
relaxation interaction between Tm and Al. We studied 


the relaxation dynamics both at and away from the res¬ 
onance using stimulated echoes. 

Away from resonance at a magnetic field of 4.56 mT, 
the echo decay was well described with a single hyperfine 
decay rate with T\ = 230 ms. At a resonance, the in¬ 
homogeneous broadening of the Tm hyperfine transition 
means that the ions have a range of different distances 
from resonance and therefore different decay rates. We 
showed that a model taking into account these different 
rates is a good fit to the observed echo decay curve for a 
minimum hyperfine lifetime of 0.5 ± 0.3 ms and a cross¬ 
relaxation interaction width of 12 ± 3 kHz. 

We also studied the spectral holeburning structure 
near resonance. For most of the resonances between Al 
and Tm, the antiholes become inverted at resonance, a 
clear sign that the hyperfine lifetime is much smaller than 
the metastable state lifetime. We modeled the structure 
of the antiholes and the central hole near resonance with 
a model that included the effect of the hyperfine inho¬ 
mogeneous broadening as well as the Zeeman structure 
of Al. In a future theoretical work, we shall model the 
cross-relaxation process, considering both the transfer of 
excitation to the neighboring Al ions and the dissipa¬ 
tion of that excitation to the environment, to see if the 
observed behaviour is consistent with the dipole-dipole 
interaction commonly invoked for cross-relaxation inter¬ 
actions between nuclear spins. 

This experimental investigation has dealt with the 
spontaneous flip flop between Tm and Al. A controlled 
flip flop can be achieved using the cross relaxation tech¬ 
nique proposed by Hartmann and HahrP. This kind of 
approach, combined with optical detection, could be ap¬ 
plied to the present system. 
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Appendix A: Solution to the 4-level rate equations 


The 4-level rate equations in Equations <[2j) are 
straightforward to solve. The solution for the popula¬ 
tion vector ft = \n\, « 2 , n- 3 , n*] 

n = c\X\ + c 2 X 2 e^ Rgt + c 3 x 3 e~ Rmt + C 4 X 4 e~ Rot (Al) 

The eigenvectors x are: 


d?i 


' 1 ' 
1 
0 

. 0 . 


(A2) 
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X 2 = 


-1 

1 

0 

0 


(A3) 


The constants c are dependent on the initial popula¬ 
tions ni(0), 712 ( 0 ), 713(0), 774(0): 


x 3 = 


Rg—Rm 

PmRm-Rg/2 

Rg—Rm 

0 

1 


(A4) 


and, for £ 4 , 


£ 4 ( 1 ) = 


(1 Pm ) Rn 

Rg - Ro 


1 Rn 


R 0 {R g - R 0 ) 

Ra 


Po R m - Ro - (1 - Pg)R, 


X 4(2) = 


®4(3) = 


1 Po Ro 

PmR m 


Rg - R 0 


flg TO 
2 1 L n 


Rg - Ro Ro(Rg ~ Ro) 


Po Rm - Rp ^ - PgRc 
1 Po Ro Rg Ro 

Rm Ro 


(1 ~Po)Ro 


£4(4) = 1 


(AS) 


(A6) 


(A7) 

(AS) 


Cl = 


?7l(0) + 77 2 (0) + 713(0) + ?74(0) 1 


C 2 = 7 : (772(0) - 774 ( 0 )) + 774(0) [ ^ ~ Pm 


Rn 


Rg Rn 


773(0) 


Po 


1 


Pg 


Ro 


Rg - Rc 


(! Po) ( 7^ Pm 


C 3 = ?7 4 (0) - 77 3 (0)(1 - Po ) 


Rn 


Rm Ro 


c 4 = n 3 (0)(l-p o ) 


Ro 


Rm Ro 


Rn 


Ro 


Rq Rm Rg Ro 


(A9) 

(A10) 


(All) 

(A12) 


For a stimulated echo, the relevant initial popula¬ 
tions 77; (0) are the populations after the first two pulses. 
Half the population resides in the non-resonant hyper- 
fine state, 77i(0) = 0.5, while the other half is divided 
between the second hyperfine state and the excited state, 
772(0) = 0.5 — do , 77.3(0) = doi where do is a measure of 
the depth of the grating. The constants c can, therefore, 


be slightly simplified: 


Cl = 2 


c 2 = do ( +Po ( o -Pg 

1 
2 

Rn 


-(! -Po) ( ^ ~P™ 


C 3 = ~do(l - Po ) 
C 4 = C?o(l ~Po) 


Rm Ro 

Ro 

Rm Ro 


Ro 


Rg ~ Ro 


Rn 


(A13) 

(A14) 


R 0 


Rg Rm Rg Rc 


(A15) 

(A16) 
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Appendix B: Structure of the A g — A e antihole 


In Section IV B| we calculated the structure of the anti¬ 
hole at A g due to the effects of the inhomogeneous broad¬ 
ening G', the resonance lineshape k a n and the hole line- 
shape C (Equation (p~5])) . This antihole was chosen be¬ 
cause it gives a simple expression, but in the experiment 
it is most often the A s — A e antihole that is seen. Be¬ 
cause the shift of this antihole is not equal to the ground 
state splitting, the equation is more complicated. 

To derive an expression for the antihole at A g — A e , we 
first make the assumption that the ground and excited 
state splittings are completely correlated so that an ion 
with a ground state splitting v c will appear in the antihole 


at a frequency v c e where e = (1 — ^). This assumption 
can be justified by looking at the holeburning structure 
in Figure |TJ the structure of the antihole at A g — A e is 
identical to that at A g , indicating that the ground and 
excited state splittings are well correlated. The lineshape 
for the antihole at A s — A e is then 


H ah (fJ c ) 



Sis + a — v A i 


x G(5i i + a — A g e)((a)da 


(Bl) 
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